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Molecular Outflows and a Mid-Infrared Census of the Massive 
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ABSTRACT 

We have observed the central region of the IR-dark cloud filament associated 
with IRAS 18507+0121 at millimeter wavelengths in CO(J=1-0), 13 CO(J=1-0), 
and C 18 0(J=1— 0) line emission and with the Spitzer Space Telescope at mid- 
IR wavelengths. Five massive outflows from two cloud cores were discovered. 
Three outflows are centered on or near an Ultracompact Hll (UCHll) region 
(G34.4+0.23) while the remaining two outflows originate from the millimeter core 
G34.4+0.23 MM. Modeling of the SEDs of the mid-IR sources identified 31 young 
stellar objects in the filament with a combined stellar mass of ~ 127 ± 27 M . 
An additional 22 sources were identified as probable cluster members based on 
the presence of strong 24 fxm emission. The total star formation efficiency in the 
G34.4 cloud filament is estimated to be ~ 7% while the massive and intermediate 
mass star formation efficiency in the entire cloud filament is estimated to be 
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roughly 2%. A comparison of the gravitational binding energy with the outflow 
kinetic energy suggests that the compact core containing G34.4+0.23 MM is being 
destroyed by its molecular outflows whereas the outflows associated with more 
massive core surrounding the G34.4 UC Hll region are not likely to totally disrupt 
the cloud. Additionally, a qualitative evaluation of the region appears to suggest 
that stars in this region may have formed in two stages: first lower mass stars 
formed and then, a few Myrs later, the more massive stars began to form. 

Subject headings: stars: formation - nebulae: HII regions - ISM: molecules 

1. INTRODUCTION 

The massive star forming region associated with IRAS 18507+0121 (hereafter IRAS 18507) 
is located in an infrared dark cloud that was mapped by Rathborne et al. (2005) in millime- 
ter and submillimeter continuum and compared to Spitzer mid-IR emission between 3.5 and 
8/xm. IRAS 18507 is roughly 11' north of the more famous ultra-compact (UC) Hll region 
complex G34.26+0.15 (Wood & Churchwell 1989; Molinari et al. 1996; Carral & Welch 
1992). The location of the infrared dark cloud relative to G34.26+0.15 is shown in Figure 1. 

Near IRAS 18507, Miralles et al. (1994) discovered a UC Hll region (G34.4+0.23) 
embedded in a 1000 M molecular cloud traced by NH 3 emission. The NH 3 emission is 
elongated in the N-S direction with a total extent of about 7', however the 1.5' resolution 
of the observations was not adequate to discern the structure of the core (Miralles et 
al. 1994). Shepherd et al. (2004: hereafter SNB04) observed IRAS 18507 in 3 mm 
continuum emission, H 13 CO + and SiO line emission and in the near-infrared. They detected 
a dense molecular core associated with the UCHll region G34.4+0.23 as well as a dense 
core and embedded protostellar object detected only in millimeter continuum emission, 
G34.4+0.23MM (hereafter G34.4MM). G34.4MM appears to be a B2 protostar surrounded 
by ~ 150 solar masses of warm gas and dust. Based on the non-detection of G34.4MM at 
near-IR wavelengths, SNB04 suggested the source may be in a very early stage of evolution. 

Rathborne et al. (2005) identified four clumps in 1.2 mm continuum emission which 
they named MM1-4 (see their figure 1). MM1 corresponds to G34.4MM, MM2 is coincident 
with the UCHll region G34.4+0.23, MM3 is about 3' north of G34.4MM while MM4 is 
about 30" south of the UC Hll region. 

IRAS 18507 was detected in a CS(2-1) survey of IRAS point sources with far-infrared 
colors suggestive of UCHll regions (Bronfman et al. 1996). The source was selected for 
further high resolution studies because of its broad line wings, a signature of current star 
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formation. By modeling HCO + , H 13 CO + , CS and C 34 S spectra obtained at an angular 
resolution of ~ 16" Ramesh et al. (1997) demonstrated that the observed line profiles can 
be explained by a collapsing hot core of about 800 M which is hidden behind a cold (~ 4 K) 
and dense (3 x 10 4 cm~ 3 ) envelope of about 200 M Q . 

The kinematic distance to IRAS 18507 has been reported to be between 3.7 and 3.9 kpc 
based on molecular line velocities that vary from 55 to 59kms _1 depending on the line 
observed and the position within the cloud (3.9 kpc: Molinari et al. 1996, Shepherd et al. 
2004, Zhang et al. 2005; 3.8 kpc: Carral & Welch 1992, Miralles et al. 1994, Ramesh et al. 
1997; & 3.7 kpc: Rathborne et al. 2005). Kinematic distances can easily have an error of 
more than 5 or even 10% due to observed deviations from circular symmetry and circular 
motion in our galaxy. Thus, a difference of 0.2 kpc is easily within the errors for this type of 
distance determination. Here we adopt the distance of 3.9 kpc to be consistent with Molinari 
et al. (1996, 1998), SNB04 and Zhang et al. (2005). If the closer distance of 3.7kpc is more 
appropriate, then mass estimates derived in this work should be decreased by about 10% 
since the mass is proportional to D 2 . 

Faiindez et al (2004) mapped the IR-dark cloud with the SEST telescope in 1.2 mm 
continuum emission and derived mass of 2000 M in the main condensation based on a dust 
temperature of 28°K. The total mass of the filament was estimated to be about 3700 M 
where 1000 M is associated with G34.4MM and 2700 M is associated with the G34.4 
UCHll region (Bronfman, personal communication). Rathborne et al (2005) mapped the 
dark cloud with the IRAM30 telescope in 1.2 mm continuum emission and derived a cloud 
mass of ~ 7500 M based on a dust temperature of 34°K, two times more mass than derived 
by Faiindez et al. In this work, we will assume the temperature and mass estimates of 
Rathborne et al. since a dust temperature of 34°K would appear to be more representative 
of massive star forming regions. However, a clear reason for why the derived temperature 
and, hence, mass estimates differ by such a large amount is not clear. 

The IRAS 18507 region is also associated with variable H 2 0, CH 3 OH and OH maser 
emission (Scalise et al. 1989; Palla et al. 1991; Schutte et al. 1993; Miralles et al. 1994; 
Szymczak et al. 2000; Wang et al. 2006; Edris et al. 2007). Molinari et al. (1996, 1998) 
observed IRAS 18507 (labeled Mol74 in their papers) and estimated a deconvolved size of 
the UCHll region of 0.7" (0.013 pc at D=3.9kpc). 

Based on observations with the NRAO 12 m telescope, Zhang et al. (2005) report that 
IRAS 18507 is not associated with molecular outflows. The lack of outflow activity contrasts 
with the presence of multiple star formation tracers (masers, dense gas, HII regions and 
millimeter continuum). 
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In this work, we present observations of the dense molecular gas near IRAS 18507 at 
millimeter wavelengths in CO(J=l— 0), 13 CO(J=1-0), and C 18 O(J=l-0) line emission to 
map the dense gas morphology and determine if there are massive outflows associated with 
IRAS 18507. Archive data from the Spitzer Space Telescope GLIMPSE survey at mid-IR 
wavelengths were studied to obtain a census of the massive and intermediate-mass star 
cluster properties over the entire infrared-dark cloud filament. Section 2 gives an overview 
of the observations, section 3 presents the results and section 4 discusses implications. 



2. OBSERVATIONS 

2.1. Observations in the 3 mm band 

Simultaneous observations in CO(J=1-0), 13 CO(J=1-0) and C 18 O(J=l-0) lines were 
made with the Owens Valley Radio Observatory (OVRO) array of six 10.4 m telescopes on 
2004 March 27 & 29 and 2004 April 16 & 18. Projected baselines ranging from about 12 to 
115 meters provided sensitivity to structures up to about 20" with ~ 5" resolution. Two fields 
were observed as a Nyquist-sampled mosaic. The observations were alternated between fields 
every 9 minutes to assure that the uv coverage in each field was similar. The total integration 
time on each field was approximately 8.8 hours. Cryogenically cooled SIS receivers operating 
at 4 K produced typical single sideband system temperatures of about 400 K. The gain 
calibrator was the quasar 1751 + 096 and the passband calibrators were 3C345, 3C 454.3 
and 3C273. Observations of Uranus provided the flux density calibration scale with an 
estimated uncertainty of ~ 15%. Calibration was carried out using the Caltech MMA data 
reduction package (Scoville et al. 1993). Images were produced using the MIRIAD software 
package (Sault, Teuben, & Wright 1995) and deconvolved with a Maximum Entropy (MEM) 
algorithm. Data analysis was performed using both the MIRIAD package and the CASA 
(Common Astronomy Software Applications) package. 

The spectral band pass for each line was centered on a velocity of 57.0 km s -1 . The 
CO images have a synthesized beam of 3.83" x 3.45" (FWHM) at P.A. -59.7°, spectral 
resolution 1.3 km s" 1 and RMS noise 60mJybeam _1 . The 13 CO images have a synthesized 
beam of 4.04" x 3.83" (FWHM) at P.A. -49.9°, spectral resolution 1.36 km s -1 and RMS 
noise 44.0 mJy beam -1 . The C 18 images have a synthesized beam of 4.04" x 3.83" (FWHM) 
at P.A. —50.1°, spectral resolution 0.683 kms -1 and RMS noise 56.0 mJy beam -1 . 
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2.2. Spitzer images 

Mid-infrared images at 3.6, 4.6, 5.8 and 8.0 fim were obtained using the Infrared 
Array Camera (IRAC; Fazio et al. 2004) of the Spitzer Space Telescope as part of the 
Galactic Legacy Infrared Mid-Plane Survey Extraordinaire (GLIMPSE; Benjamin et al. 
2003). Image mosaics were created by the GLIMPSE pipeline after image artifacts such 
as cosmic rays, stray light, column pull-down, and banding were removed]]] Point sources 
were extracted from the GLIMPSE images using a modified version of Daophot (Stetson 
1987). The astrometric uncertainty for point sources is ~ 0.3". The GLIMPSE 5<r point- 
source sensitivities are approximately 0.3, 0.3, 1.2, and 0.7 mJy in the four IRAC band 
passes, respectively. Point source sensitivities were generally higher in areas of bright diffuse 
background emission. The fluxes of point sources that are detected at > 5a twice in one 
IRAC band or once in each of two bands are entered in the GLIMPSE Point Source Archive. 
Each Archive entry includes the corresponding J, H, and K s fluxes from the 2MASS Point 
Source Catalog (Skrutskie et al. 2006) if that source was detected in 2MASS. 

The Multi-band Imaging Photometer for Spitzer (MIPS; Rieke et al. 2004) was used to 
image the entire GLIMPSE survey area at 24 and 70 /xm as part of the MIPSGAL survey 
(Carey et al. 2005). Fluxes at 24 fim for point sources within the G34.4 IR dark cloud 
were extracted from the frame using the GLIMPSE point-source extractor. The relative 
uncertainty between two sources in a MIPSGAL image is roughly 0.3 to 0.5 pixels while the 
absolute positional uncertainty is ~ 1". 



3. RESULTS 

Figure 1 presents a four-color (3.6, 4.6, 5.8, & 8.0 jttm) GLIMPSE image of the G34.4 
environment in Galactic coordinates (North points to the upper left corner). The "cometary" 
UCHll region G34.26+0.15 is located on the right side of the image. An IR-bright nebula 
extends below the G34.26 UC Hll region. The UC Hll region contains a cluster of OB 
stars, the most luminous of which is an 06.5 spectral type (Wood & Churchwell 1989). The 
nebula is bright in the IRAC bands primarily due to PAH emission features stimulated by 
UV radiation from the cluster. The subject of this paper is the dark filament in the upper left 
which is forming early B stars. The mid-IR emission from the UCHll region G34.4+0.23 is 
significantly less than that of G34.26, as expected due to the lower luminosity of the central 
source. Figure 2 presents a close-up of the IR-dark filament, now displayed in Equatorial 



1 Details of the data processing can be found at 



http: / / www.astro.wisc.edu/glimpse/docs.html 
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coordinates with North pointing upwards. Along the dark filament, North and South of 
G34.4, are embedded sources seen predominantly by their 4.5 /zm (green) emission. Our 
OVRO CO mosaic is centered on the southern core near the UCHll region G34.4 and the 
millimeter core G34.4 MM located ~ 40" north of the UC Hll region. 

Figure 3 presents the red and blue-shifted emission near G34.4MM and G34.4 along 
with a 2 /im K band image from SNB04. Neither the millimeter core nor the UC Hll region 
is visible in the K band image. Both regions show complex high-velocity structure suggesting 
the presence of multiple outflows. 

Figure 4 shows the integrated emission (moment 0) in C 18 and 13 CO along with the 
red and blue-shifted 13 CO emission. The bottom two panels in figure 4 show moment 1 
maps of C 18 & 13 CO in which the brightness of the color is proportional to the intensity 
in the moment image. This type of image provides some insight into how the dense cloud 
structure and kinematics are being affected by the molecular outflows. 

Figures 5, 6 & 7 present channels maps of CO, 13 CO and C 18 emission, respectively. 

3.1. Molecular outflows 

At least two outflows can be seen emanating from G34.4MM (labeled A and B in the 
upper left panel of Fig. 3) while at least three outflows originate from the southern core near 
the G34.4 UCHll region (labeled C, D and E). 

The velocity structure of the outflows from the G34.4MM region is consistent with a 
single outflow dominating the red-shifted kinematics (outflow G34.4:A). There is a weak 
blue-shifted counter flow extending to the NE that can be traced to G34.4MM. The red- 
shifted lobe of outflow G34.4:A appears to have little or no associated blue-shifted emission 
suggesting that it is oriented outside of the plane of the sky (the blue-shifted emission in the 
red lobe appears to be from the collimated flow G34.4:E coming from the southern region). 
The full-opening angle of the red-shifted lobe is ~ 43°. 

Marginal evidence for an outflow perpendicular to outflow G34.4:A is seen in the red 
and blue-shifted emission clumps along the axis labeled 'B'. No source has been detected at 
millimeter wavelengths which could drive outflow G34.4:B. The source could be embedded 
in G34.4MM and not resolved at the resolution of ~ 5" of SNB04. 

The southern molecular core containing the G34.4 UC Hll region contains three over- 
lapping flows seen in the CO emission (flows 'C, 'D' and 'E'). Blue-shifted CO, 13 CO and 
C 18 gas is seen toward the G34.4 UC Hll region extending to the West (flow G34.4:C). This 
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outflow is seen most clearly in 13 CO and C 18 emission in Figure 4. 13 CO blue-shifted emis- 
sion is traced back to the Hn region (top right panel in Fig. 4) while the bottom two panels 
illustrating the velocity structure in the dense gas show that blue-shifted emission extends 
from the UCHll region to the west. A possible red-shifted outflow lobe east of the UC Hn 
may be visible in 13 CO emission in the top right panel in Fig. 4 although the identification 
is not definitive. 

The NE-SW flow (G34.4:D) appears to be centered about 10" SE of the UC Hn region, 
just off the peak of a dense C 18 clump. The flow appears to be relatively symmetric in 
the red and blue-shifted gas distribution although the actual morphology is considerably 
confused due to overlapping outflows. If the flow is symmetric, then one would expect the 
driving source to be located near the small, filled triangle at 18 h 53 m 19.0 s , 01 o 24'35" along 
the G34.4:D outflow axis in Figs. 3 & 4. 

The SE-NW flow (G34.4:E) is seen clearly in both CO and 13 CO emission. If the driving 
source is located between the red- and blue-shifted emission along the outflow axis, it would 
be centered about 5" S-SW of the UCHll region (indicated by a small, filled triangle at 
18 h 53 m 18.6 s , 01°24'43" along the G34.4:E outflow axis in Figs. 3 & 4). The outflow is not 
symmetric: the red-shifted emission lobe extends about 20" to the SE while the blue-shifted 
lobe extends roughly 60" to the NW, overlapping the blue-shifted lobe from the G34.4 UC Hn 
region flow and the red-shifted lobe of flow G34.4:A farther north. 

Outflow G34.4:A has a position angle of +58° (measured from North to the axis of 
the blue-shifted lobe) while G34.4:B P.A. = -38° and the difference in orientation between 
outflow axes is 96°. The positions angles of the three outflows near the G34.4 UC Hn region 
are: G34.4:C P.A. = -100°; G34.4:D P.A. = +38°; G34.4:E P.A. = -34° and the maximum 
difference in orientation between the axes is 138°. 



3.2. Mass &: Kinematics of the Outflows 

The mass associated with CO line emission is calculated following Scoville et al. (1986). 
The CO excitation temperature is taken to be 34 K which is the dust temperature in MM1 
derived from fits to the SED (Rathborne et al. 2005). We assume the gas is in LTE with 
[CO]/[H 2 ] = 10~ 4 , [CO]/[ 13 CO] = 51 and [CO]/[C 18 0] = 378 at the galacto-centric distance 
of 5.8 kpc (Wilson & Rood 1994). 

The CO optical depth at red and blue-shifted velocities is estimated from the CO/ 13 CO 
ratio of spectra taken at four positions in the mosaic image that have been convolved with a 
10" beam. The open triangles in the top, center panel of Fig. 4 show the locations where the 
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spectra were taken, while Figure 8 shows the spectra taken at those positions. In channels 
where no 13 CO is detected, we assume that the CO is optically thin. The CO and 13 CO line 
profiles show a dip near ^lsr indicating that the lines suffer from self-absorption and/or the 
interferometer is missing extended flux at lower velocities. 

Because multiple, overlapping flows are present, it is not possible to obtain a reasonable 
estimate of the inclination of each flow. Thus, we assume an inclination angle of 45° which 
minimizes errors introduced by inclination effects for flows that are not near the plane of the 
sky or along the line-of-sight. Table 1 summarizes the physical properties of the molecular gas 
in the combined outflows originating from G34.4MM and those centered near the UCHll 
region. The total flow mass Mf is given by ^ where Mj is the flow mass in velocity 
channel i corrected for optical depth. The momentum P is given by Y MiVi and the kinetic 
energy E by \ Y Mivf where Vi is the central velocity of the channel relative to v lsr- 

The dynamical timescale t d is calculated using Rf/ < V >, where the intensity- weighted 
velocity < V > is given by Pj (Y (Cabrit & Bertout 1990) and Rf is the average outflow 
radius of the flows. This method of calculating t d is really only correct for a single outflow; 
however, an independent outflow age cannot be determined for both flows in the cores since 
they cannot be separated near the central position. Because of this uncertainty, the flow 
dynamical timescale could easily be in error by a factor of 2 to 5. 

The mass outflow rate Mf is ^ Mi/t d , the force F is P/t d and the mechanical luminosity 
L is E/td- Since these values are dependent on t d x uncertainties in t d are directly applicable 
to Mf, F, and L. 

The combined outflows from the G34.4 core have a total mass of 111 M and energy, E 
= 7.8 x 10 46 ergs. The G34.4MM outflows have a combined mass 34.8 M , E = 4.9 x 10 46 ergs. 
Despite the uncertainties, the flow masses and energies are consistent with those for outflows 
driven by young, early B stars. 

3.3. The dense gas cores 

The C 18 and 13 CO moment images (Fig. 4) show clearly that the two cores have 
very different characteristics. The southern core centered on the G34.4 UCHll region is 
extended, with a largest spatial scale of about 1.4 pc at a distance of 3.9 kpc. The northern 
core, centered on G34.4 MM, is significantly more compact with higher peak emission than 
the southern cloud. The diameter of the G34.4MM core is about 0.2 pc and is surrounded 
by faint, diffuse 13 CO emission that extends out to about lpc in size. 
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The total mass of the 2 cores as measured from C 18 is given in Table 2. C 18 is 
assumed to be optically thin although this is not likely to be true near vlsr & t the peak of 
the line profile. The extended, southern cloud in which the UCHll region is embedded is 
almost 10 times more massive (692 M ) than the G34.4MM core (75 M ). 

As Fig. 8 illustrates, the C 18 line appears to be suffering from missing flux due to 
the lack of short baselines. This is most apparent in the spectra taken SE of the UC Hll 
region position where C 18 is more extended and there is a clear double peaked line profile. 
The C 18 line is single peaked closer to the cloud core suggesting that the interferometer is 
recovering most of the flux density and we are not seeing significant self-absorption. Thus, 
the mass estimates based on C 18 are lower limits. 

The G34.4 and G34.4MM molecular cores have a projected separation of about 40" 
(0.75 pc at 3.9 kpc). Their LSR velocities as measured by the peak of the C 18 line differs 
by 1.4 km s -1 . 

3.4. Mid-Infrared Sources: GLIMPSE & MIPSGAL 

G34.4MM and the G34.4 UCHll region are embedded in an extended, filamentary 
IR-dark cloud complex seen against the diffuse PAH background emission which is visible 
primarily in the Spitzer 8.0 /zm band (Figures 1 & 2). The G34.4MM core is associated with 
diffuse emission in all 4 IRAC bands, but is particularly noticeable at 4.5 /im. The [4.5] 
band contains the H I Bra line along with some molecular forbidden lines, notably the CO 
1-0 P(8) fundamental band-head, but, unlike the other 3 bands, does not contain any PAH 
emission feature. The [4.5] band also contains several H 2 lines (e.g. the H 2 0-0 S(9) line 
at 4.694 fim) which trace shocks in some massive outflows (e.g. DR21: Smith et al. 2006; 
W75N: Davis et al. 2006). Thus, the [4.5] band appears to readily trace ionized and/or 
shocked gas. 

Figure 9 shows the IRAC 4.5 /im image in greyscale with 13 CO contours overlayed for 
reference. The plus symbols represent the locations of the UCHll region and G34.4MM 
core. The open diamonds just south and south-east of the UCHll region represent the 
predicted locations of sources driving outflows G34.4:E (18 h 53 m 18.6 s , 01°24'43") and G34.4:D 
(18 h 53 m 19.0 s , 01°24'35") based on CO and 13 CO emission. While the predicted position of 
the YSO that was identified based on the CO outflow morphology should be considered no 
better than an educated guess, it is interesting that each position is associated with a mid-IR 
source to within a few arcsec. The characteristics of these sources along with other sources 
in the cluster are discussed below. 
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The spectral energy distributions (SEDs) were analyzed by fitting model SEDs to each 
observed source from a large pre-computed grid of YSO models (Robitaille et al. 2006, 2007). 
The grid consists of 20,000 2-D YSO radiation transfer models spanning a complete range of 
stellar mass and evolutionary stage. Each YSO model outputs SEDs at ten viewing angles 
(inclinations), so the fitter actually has 200,000 SEDs to choose from. The YSO models 
also output SEDs for each of 50 aperture sizes. Thus, if the distance to the star formation 
region is known, only the model SEDs with the appropriate aperture size compared to the 
observed photometry are used. A grid of stellar atmospheres were also incorporated into the 
fitting algorithm in order to facilitate the separation of main-sequence and evolved stars from 
YSOs (Brott & Hauschildt 2005; Kurucz 1993). The fitting algorithm includes a foreground 
extinction component to the model which allows one to easily distinguish highly reddened 
stars from YSOs which have an additional emission component in the mid-infrared. 

The 24 /im fluxes were cross correlated with positions of the 2199 GLIMPSE Archive 
sources lying within a 3.9' radius of the G34.4 UCHll region to account for observational 
errors in addition to photon counting statistics. For the purposes of the fitting, all flux 
uncertainties that fell below the 10% level were reset to a minimum value of 10%. The SED 
fitter was then run on this combined 8-band source list, fitting only the fluxes of sources 
detected in 4 or more bands. Sources with SEDs well fit by a stellar atmosphere are not 
likely to be cluster members. Sources with a x 2 per flux data point (xV^data) of 4.8 provided 
a reasonable dividing line between sources that, by eye, appeared to be well- fit by stellar 
atmospheres and those that were not well-fit due to IR excesses. Of the total 2199 sources 
in the GLIMPSE Archive in this region, 868 were detected in four or more bands. Of those, 
829 sources were well-fit by stellar atmospheres and thus are likely foreground or background 
stars rather than YSO cluster members. The remaining sources were fit with YSO models 
and 27 of these met the xV-^data < 4.8 criterion. Of the 9 sources that were not well-fit either 
by stellar atmospheres or YSO SEDs, 1 appears to be a YSO and 8 are probably Archive 
sources with improperly extracted fluxes in one or more bands or there was a mismatch of 
GLIMPSE with 2MASS or MIPS sources. 

An overview of the results of the SED fitting procedure is presented in Figure 10 and 
Tables 3 & 4. Most of the filamentary dark cloud is visible in the three-color image (Fig. 10) 
in which red corresponds to MIPS 24 /im, green to GLIMPSE 8 /im and blue to GLIMPSE 
4.5 /im. The G34.4 complex is located near the center of Fig. 10 (in the region outlined 
with a large, white circle). Point sources are marked by small circles color-coded by source 
type as identified by the Fit flag in Table 3. Sources that were well-fit as YSOs are marked 
by green circles and given a Fit flag of 1 in Table 3. The yellow circle identifies a probable 
YSO that was not well- fit (Fit = 2). Blue circles identify sources that were well-fit by 
stellar photospheres. MIPS 24 /xm point sources that are detected in fewer than 3 of the 
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2MASS+IRAC bands and hence were not run through the fitter (Fit=4) are identified by red 
circles; the majority of these are probably YSOs. The cyan circle (source 32 located between 
sources 8 and 10 in Fig. 10) marks the Archive source that coincides with the proposed origin 
of outflow G34.4:D, it is detected only at 3.6 and 4.5 /zm, and hence was not run through the 
fitter (Fit=3). The area covered by the OVRO observations is contained primarily within 
the large white circle, which is referred to as "region C" in Table 3. "Region N" is identified 
by the box north of region C, while "region S" is identified by the box to the south. The 
G34.4 UC HII region appears white in this image because it is bright in diffuse emission at 
all wavelengths. There are two prominent MIPS sources: G34.4MM and another, fainter 
source located about 2' north of G34.4MM at position 18 h 53 m 20.602 s , +01°28 / 25.61". 

Table 3 presents a summary of the positions and fluxes in each band for all identified 
probable YSO cluster members in and near the dark cloud. Table 4 presents the best fit model 
parameters for each YSO with enough flux points to be fitted. Output model parameters 
are: interstellar extinction to the source, Ay; stellar mass, M*; total luminosity of the YSO, 
L+; envelope accretion rate, M env ; and inclination, i. Note that Ay does not include the 
extinction produced by a circumstellar disk or envelope and thus it represents a lower limit 
for YSOs with significant circumstellar material. In general, each source can be well fit by 
multiple YSO SED models and this generates a corresponding range of the best-fit model 
parameters. For well-fit YSOs, results of all fits with Ax 2 /iV data < 1 relative to the best- 
fit model are reported. E.g., for a typical source detected in iV data = 5 bands, all fits with 
Ax 2 > 5 with a formal probability relative to the best fit of less than exp(— A% 2 /2) = 0.08 are 
excluded. The cumulative probability distribution for each parameter is then constructed 
using the Gaussian- weighted probability of each fit relative to the best fit. A range of 
parameter values that has a greater than 95% probability of containing the actual value for 
the source is then determined and reported in Table 4. 

The central region of the filament: 

In the region for which CO outflow data is determined, 7 YSOs with good fits have been 
identified (sources 6-12), along with one poorly-fit, possibly massive YSO (source 31 in tables 
3 & 4) located on the northeast edge of the G34.4 MM 24 jum emission. Source 12, the highly 
luminous, the mid-IR counterpart to G34.4 MM, is most likely the driving source of outflow 
G34.4:A. Marked by a larger green circle in Figure 10, source 12 is partially resolved in the 
GLIMPSE images (this is seen most clearly in the 4.5//m greyscale image in Fig. 9). The 
morphology of source 12/G34.4MM is consistent with two dust emission lobes separated by 
an absorption lane due to a disk where the lobes are in rough alignment with the axis of 
outflow G34.4:A. 

The fluxes in the two lobes of source 12/G34.4MM in each GLIMPSE band was summed 
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and combined with the point-source fluxes from MSX 21.3 /iin (Price et al. 2001), MIPS 24 
/im and MIPS 70/im to model the source SED across the mid-IR. The 10 best fits to source 
12/G34.4MM are plotted in Figure 11. The fit parameters are well-constrained and the best 
fit has x 2 = 6-62 which is significantly lower than the second-best fit with y 2 = 20.79. The 
best fit parameters suggest a massive YSO with M* ~14M , ~ 9,400 L . The source is 
deeply embedded, behind greater than Ay ~ 40 mag of extinction and it is young with a 
high envelope accretion rate of M CDV ~ 10~ 3 M yr _1 . 

The inclination angle of the best-fit model is i = 56°, consistent with the orientation 
of outflow G34.4:A which has a wide opening angle of 43° but does not intersect the plane 
of the sky. This interpretation assumes that a single YSO dominates the flux density and 
geometry of the mid-IR emission. Given that two outflows have been identified emanating 
from the G34.4 MM core and the far-IR fluxes exceed those predicted by a single YSO SED, 
more than one YSO is likely present. However, the resolution of the GLIMPSE images is 
not adequate to clearly distinguish between sources that are likely to drive outflows G34.4:A 
and G34.4:B. 

For the remaining six well-fit YSOs in region C, the best fit models are presented in 
Figure 12a while the distribution of their corresponding best-fit masses and luminosities are 
shown in Figure 12b. Of these, two YSOs (sources 10 & 11) are nearly coincident with the 
UC Hll region. 

Source 11 is a luminous YSO located near the projected center of the UC Hll region and 
appears to be the most likely candidate to drive outflow G34.4:C. The GLIMPSE position 
differs from the 6 cm ionized gas peak by ~ 1.3". This could be due to differences in 
astrometry between the radio and IR data; however, the discrepancy exceeds the GLIMPSE 
astrometric uncertainty of 0.3" and the Very Large Array estimated uncertainty of ~ 0.1". 
It is more likely that the centimeter peak tracing the ionized gas is offset from the actual 
YSO. This could easily happen if the ionized gas morphology closely follows the confines of a 
bipolar outflow cavity. It is also possible that the center of the source is not detected in the 
IRAC bands due to high extinction or the bright diffuse background of the UCHlI region 
itself. We may only see the warm outflow cavity associated with the blue-shifted lobe in the 
mid-IR that is displaced from the centimeter continuum emission. 

Both sources 10 and 11 suffer from confusion with the UCHll region at IR wavelengths 
longward of the IRAC bands. The UC Hll region emits very strongly in all 4 MSX bands 
and in MIPS 24 and 70 /jm. These fluxes were used as upper limits for the fits to sources 
10 and 11, but they are very high limits (see Figure 12a), and the fits to both sources are 
poorly constrained. The best fit results for source 11 suggest M± =2-10 M and L* =200- 
5,200 L , but the probability distribution of fit parameters shown in Figure 12b is skewed 
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and allows for higher masses and luminosities up to ~ 20 M , L+ ~ 45, 000 L Q (albeit at 
a formal probability level of only a few percent). SNB04 determined a spectral type of BO. 5 
(L* ~ 40, 000 L Q ) for the central star of the G34.4 UC Hll region assuming a single zero-age 
main sequence (ZAMS) star is producing the observed Lyman continuum flux. Thus, if 
source 10 is the ionizing star of the UC Hll region, the higher values of mass and luminosity 
are more plausible. 

Due to its physical location between the red and blue-shifted lobes of outflow G34.4:E, 
source 10 appears to be the most probable driver of this outflow. The best fit results suggest 
that it is a YSO with M* =0.5-10 M and L* =45-4,500 L . Again, the association of this 
YSO with strong outflow activity suggests that the upper end of this range of mass and 
luminosity may be more appropriate. 

Source 32 (cyan circle), the potential driver of outflow G34.4:D, is detected in only 
two IRAC bands ([3.6] and [4.5]) at the 5a level sufficient for inclusion into the GLIMPSE 
Archive. Based on these two bands, source 32 has a red color of [3.6] — [4.5] = 1.3 which 
is consistent with it being a YSO (Robitaille et al. 2006). The lack of GLIMPSE Archive 
fluxes in bands [5.8] and [8.0] means that this source was not detected at the 3a level in 
those two bands, corresponding to magnitudes fainter than 13.0 in each bano^J The colors 
of YSO models for very young sources predict [3.6] — [5.8] > 0.4 and [5.8] — [8.0] > 0.4, 
corresponding to magnitudes fainter than 13.6 and 13.2 in bands [5.8] and [8.0], respectively 
(Robitaille et al. 2006). Thus the non-detection in bands [5.8] and [8.0] is consistent with 
the interpretation that source 32 is a YSO. 

Source 32 is also not detected in the MIPS 24 fim band. However, given the proximity 
to the very bright UCHll and lower resolution of MIPS compared to IRAC, this source 
is confused with the PSF wings of the UCHll region. Based on an estimate of the MIPS 
sensitivity in a confused region, a 3a upper limit of 50 mJy (magnitude 5.4) is expected 
for the 24 /zm flux. This could potentially allow a very red [3.6] — [24] w 8.6 color for this 
source. We therefore conclude that the mid-IR colors of source 32 are consistent with those 
for a YSO, however the strongest evidence for this remains the fact that it lies along the CO 
outflow axis of G34.4:D near where one might expect the driving source to be located. 

In summary, 9 YSOs have been identified and characterized in region C near G34.4MM 
and the UC Hll region. Four or five are probably massive with M+ > 8 M while the others 
are low-to-intermediate mass with M* between 0.5 and 8M . 



2 See the GLIMPSE Quality Assurance Document at http://www.astro.wisc.edu/glimpse/GQA- 
master.pdf 
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The northern filament: 

North of the G34.4MM and UCHll region complex (in region N) 17 YSOs have been 
successfully identified based on fits to the SEDs. Most of these YSOs closely trace the 
IR-dark cloud structure. At the extreme north end of the field there is another 24 /iin bright 
source (source 29 in Tables 3 & 4) that may be similar to but less massive than G34.4 MM. 
The SED is well-fit with X 2 = 17.39 suggesting a YSO with M* ~4.8-6.8M and L* ~170- 
340 L . The remaining 16 YSOs in region N are less massive than Source 29 with mass 
estimates between about 0.5 and 6.5 M Q . 

The southern filament: 

The southern part of the IR-dark filament (region S) contains 5 sources with SEDs that are 
well-fit as YSOs. Five relatively faint 24 fj,m sources that are likely to be YSOs are also 
scattered throughout this region. Sources in this southern region are more distributed and 
are not spatially well-correlated with the dark cloud filament. This is in contrast to the 
YSOs in the northern region where they follow the cloud boundaries reasonably well. Mass 
estimates range from about 1 to 7M . 

4. Discussion 

The IR-dark cloud associated with IRAS 18507+0121 has a filamentary structure that 
spans more than 9' on the sky (> 10 pc at 3.9 kpc). Splitting the cloud into three regions 
(south, central and north) one finds that the southern region is associated with mostly low 
to intermediate-mass YSO candidates (10 in total) that are not spatially well-correlated 
with the dark cloud filament. The northern segment of the filament contains 23 low to 
intermediate-mass YSO candidates that tend to trace the IR-dark cloud structure. The 
central region of the filament is in the process of forming four or five mid- to early-B stars 
with M* > 8M . Four more sources are likely low-to-intermediate mass with between 
0.5 and 8M . Together with an additional 11 YSO candidates seen at 24 /xm this central 
region harbors a relatively compact cluster with a total of 20 YSO candidates. 

4.1. GLIMPSE-detected source characteristics 

Bright emission in the 4.5 /zm Spitzer IRAC can be a strong indicator of molecular 
outflow (e.g. Qiu et al. 2007, Smith et al. 2006, Noriega-Crespo et al. 2004). The most 
likely reason for this is the presence of H 2 emission lines and/or emission from the CO 1- 
P(8) bandhead. All four IRAC bands contain H 2 emission lines, but these features are 
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narrow and hence can only significantly enhance the brightness in the broadband images if 
the lines are extremely bright relative to the continuum. Because H 2 features are distributed 
throughout the IRAC bands, strong H2 lines will not cause a greater enhancement of the 
flux in band [4.5] than in bands [5.8] and [8.0] (see e.g., Smith et al. 2006). In contrast, 
the CO 1-0 P(8) emission bandhead is a broad feature that occupies the entire [4.5] filter 
and can become bright in the presence of shocked, molecular gas. There is not significant 
CO bandhead emission in the remaining IRAC bands. Thus, sources with excess emission 
in band [4.5] relative to bands [5.8] and [8.0] appear likely to have CO bandhead emission 
resulting from shocked, molecular gas. This does not preclude the possibility that significant 
H 2 emission may also contribute to the brightness in all IRAC bands. 

Three of the 29 YSOs listed in Table 3 that were detected in all four IRAC bands have 
noticeable excess emission at 4.5yum and strong flux at 24/iin: sources 6 & 12 in region C, and 
source 29 in region N. Because of this, SED fits used the [4.5] fluxes as upper limits only since 
the YSO SED models do not incorporate molecular line emission. These sources are likely 
to have significant shocked, molecular gas and are undergoing outflow & accretion. Source 
31, although it was not well-fit, should probably be considered in this group as well. It was 
not well-fit because it has a strong 4.5/xm excess and it was only detected in the four IRAC 
bands. The SED (not shown) looks similar to G34.4MM itself. Further, its non-detection 
at 24/iin is almost certainly due to its proximity to G34.4MM. 

The most prominent of the 4.5/im excess sources is source 12, the mid-IR counterpart to 
G34.4 MM and outflow G34.4:A. The association of this source with the CO outflow favors 
the interpretation that the 4.5/im excess is due to shocked CO emission. 

The potential drivers of outflows G34.4:C (source 11) and G34.4:E (source 10) do not 
show similar 4.5/im excess fluxes. A possible reason for this may be that they are somewhat 
more evolved than source 12. In very young sources, particularly those with M* > 2 M , 
high accretion rates create a cool, distended photosphere that does not emit strongly enough 
in the UV to excite the PAH emission features present in the IRAC [3.6], [5.8] and [8.0] 
bands. In more evolved sources, however, PAH emission can hide the contrast between [4.5] 
and the other 3 bands. Alternatively, as discussed above, both sources 10 and 11 suffer from 
confusion with the UC Hll region which emits strongly in all four bands and this could also 
mask moderate enhancements due to CO bandhead emission. 

As a final note, three sources show marginal excess at 4.5/im relative to bands [5.8] and 
[8.0]: source 7 in region C and sources 20 & 26 in region N. These sources do not appear 
to have strong emission at 24/iin, so they may be a different type of source than those with 
strong 4.5 & 24yum excesses. 
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4.2. Timescale for Massive Star Formation 

There are two cloud cores in the central region of the filament: one more massive core 
in the south surrounded by a more diffuse molecular cloud that contains the UC Hll region 
and one that is more compact containing source 12/G34.4 MM. The cores are separated by 
about 0.75 pc. Each cloud core is producing a small cluster of massive flows that look similar 
and are of similar ages although the northern outflows associated with G34.4 MM have about 
one third the mass as the southern flows. 

Due to the shorter formation timescales (^,10 5 years), massive (proto)stars begin to 
burn hydrogen even while they continue to accrete material. The generation of a UCHll 
region identifies which sources have reached the main sequence and begun hydrogen burning 
while the detection of a molecular outflow traced to the source identifies whether it is still 
accreting (and hence still powering an outflow). However, the addition of material onto the 
stellar surface implies that the source will move up the main sequence. In this sense it is 
still a protostar because it has not reached its final mass. 

The central source in the G34.4 UCHll region (source 11) has a spectral type of BO. 5 
while the spectral type of source 12/G34.4MM is probably closer to a B2 star (this work & 
SNB04). High- velocity CO emission can be traced to both sources suggesting that they are 
still accreting even though they have reached the main-sequence. 

The characteristic time scale of the outflow/infall phase from early B stars is roughly a 
few xlO 5 years (Richer et al. 2000) while the evolutionary timescale for a BO. 5 star to reach 
the main-sequence is roughly 10 5 years. Since ionized gas is seen toward both cores and the 
outflow characteristics are similar, it seems reasonable to assume that they began forming 
massive stars at approximately the same time: about 10 5 years ago. 

SNB04 detected a lower mass population of this cluster (e.g. a few solar masses or less) 
at near-infrared wavelengths (see their Fig. 6). Their detected cluster population was more 
distributed, less embedded and less massive than those detected in this work. They found 
that this sub-population is less than 3Myrs old but probably more than 0.3Myrs old. 

There are numerous problems associated with these age estimates and they could be in 
error easily by a factor of a few. However, let us take them at face value for now. SNB04 
discovered a distributed population of low to intermediate mass stars that may have formed 
a few Myrs ago. The more massive stars detected in this work are more tightly associated 
with the dense molecular gas (as expected) and appear to be somewhat younger, more like 
10 5 years old. This suggests that the star formation in this cold, dark cloud may have formed 
in two stages: first lower mass stars began to form, then, a few Myrs later, the more massive, 
early B stars began to form. 
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Such a situation is predicted when clouds have magnetic fields that are stronger than 
needed to provide support against gravitational collapse and the core mass to magnetic 
flux ratio (M/$) is sub-critical. Thus, low-mass stars would form in originally highly 
magnetically sub-critical clouds, with ambipolar diffusion leading to core formation, quasi- 
static contraction of the cores, and eventually, massive star formation (see e.g. the review 
by Ward-Thompson et al (2007) and references therein). If this scenario is correct, then the 
cold, dark filament associated with IRAS 18507 may provide an example of a cloud with 
delayed massive star formation due to an initially high magnetic field and sub-critical M/$. 



4.3. Star Formation Efficiency 

The set of best-fit SED models for YSO candidates can be used to construct a x 2 - 
weighted probability distribution of the best-fit model parameters. We sum the mass 
probability distributions for each source to construct mass functions for all well-fit YSOs 
(Fit=l) in Tables 3 & 4. The total mass of the GLIMPSE-detected YSOs in this IR-dark 
cloud is then 127±27M . 

The total stellar mass of the central (sub)cluster containing G34.4MM and the UCHll 
region is 48 ± 20 M (or ~ 38% of the total detected stellar mass). About 14 M is in source 
12/G34.4MM, 10 M is in source 31 (the yellow source near G34.4MM in Fig. 10), and the 
remaining 24 M of stellar mass resides in the central star producing the UC Hll region and 
its surrounding cluster. These rough stellar mass estimates are incomplete for M* < 3M . 
They also do not include the additional 20 YSO candidates that were detected at 24 /zm but 
in fewer than three IRAC bands. Thus, the stellar mass estimate should be considered a 
lower limit. 

Using the initial mass function (IMF) of Kroupa et al. (2001) (where the IMF slope 
follows a Salpeter IMF for M* > 1M and it becomes shallower for 0.1 < < 1M ) one 
can estimate the expected stellar mass for M* < 3M . Assuming 127 ± 27 M of mass in 
stars more massive than 3 M , roughly 413 ± 103 M of mass is expected to be in the form of 
stars less massive than 3M . Thus, the total stellar mass content of this cloud is expected 
to be ~ 540 ± 130 M . 

This stellar mass estimate is a little more than half of the derived mass of the 
ammonia hot core of 1000 M surrounding the UC Hll region (Maralles, Rodrigues & 
Scalise 1994). And it is about 7% of the cloud mass of 7500 M derived from 1.2mm 
continuum emission (Rathborne et al. 2005). The total star formation efficiency is then 
M*/(M doud + M*) = 540/(8040) ~ 0.07 or 7%. This should be considered an upper limit 



- 18 - 



because 1.2 mm continuum emission traces warm dust (typically 30-50 K) which is heated 
predominantly by massive stars - it does not effectively trace the cold gas and dust (5-10 K) 
which may represent the bulk of the total cloud mass. 

We can also obtain an estimate of the massive and intermediate mass star formation 
efficiency in the warm gas traced by 1.2 mm continuum emission and in the densest part 
of the hot core. Assuming a warm cloud mass of 7500 M and a total mass of stars with 
M± > 3M = 127 M , the massive star formation efficiency is roughly 2%. 

4.4. Can the outflows destabilize the cloud cores? 
The G34.4 cloud core surrounding the UC Hll region: 

Maralles, Rodrigues & Scalise (1994) estimated a mass of the G34.4 cloud core surrounding 
the UC Hll region to be ~ 1000 M based on NH 3 observations. They only detect emission 
near the UCHll region, not near the G34.4MM core. The C 18 mass estimate of 693 M 
compares reasonably well with that of NH3 given the uncertainties inherent with both 
estimates and the fact that our C 18 map is missing extended emission. 

If we take the mass of the G34.4 core surrounding the UCHll region to be M c i OU( i = 
1000 M we find that the core has about nine times more mass than in the high- velocity 
molecular outflows. The gravitational binding energy of the cloud, can be calculated from 
GM^ [oud /cir, where C\ is a constant which depends on the mass distribution (ci = 1 for 
p oc r~ 2 ). We find that the gravitational binding energy is 2.4 x 10 47 ergs. Roughly 11% of 
the molecular cloud core is participating in the outflow and the combined outflow energy is 
roughly one third the gravitational binding energy of the cloud. The G34.4 outflows near 
the UC Hll region are injecting a significant amount of mechanical energy into the cloud core 
and may help prevent further collapse of the cloud. However, they are not likely to totally 
disrupt the cloud. As long as cloud material remains, this leaves the door open for possible 
future episodes of star formation. 

The compact G34.4MM cloud core: 

SNB04 estimate the mass of the G34.4MM core based on 2.6 millimeter thermal dust 
emission to be between 150 and 650 M with 240 M being the value derived for an opacity 
index, (3 of 1.5 and a temperature of 50 K. Our estimated mass based on C 18 emission is 
significantly less than this (~ 75 M ) suggesting that our C 18 maps are missing extended 
emission or the assumptions used to estimate the mass associated with warm dust emission 
in the G34.4 MM core were not correct. Thus, for example, if (3 = 0.7 and T d = 70 K then the 
mass of the millimeter core would be ~ 75 M which is consistent with the mass derived from 
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C 18 0. Such a low opacity index and high temperature could be caused by the presence of 
strong shocks through most of the core mass which would prevent large grains from forming 
and maintain the bulk of the gas at higher temperatures. Indeed, this could be possible 
given the small size of the core (0.2 pc) and the fact that two outflows with a combined mass 
of nearly 35 M are blasting their way out of the core. 

If we assume a mass of 75-240 M for the G34.4MM core, as derived from the C 18 
and 2.6 mm continuum observations, then the gravitational binding energy of the cloud is 
5.6 x 10 45 ergs to 5.7 x 10 46 ergs. The kinetic energy of the outflowing gas is between 1.2 to 
nearly nine times the gravitational binding energy of the cloud core. Thus, the G34.4MM 
flows appear to have more energy relative to the core and are more likely to disrupt the 
compact core. Therefore we may be seeing the first and only episode of massive star formation 
in this compact region. 

5. Summary 

We have observed the central region of the IR-dark cloud associated with IRAS 18507+0121 
in CO J=l-0, 13 CO(J=1-0) and C 18 O(J=l-0) with the Owens Valley Radio Observatory 
at ~ 4" resolution. Five massive outflows from two cloud cores were discovered. A total of 
146 M of material is participating in the combined outflows and injecting ~ 1.3 x 10 45 ergs 
of kinetic energy. We also used archived Spitzer data from the GLIMPSE survey to gain an 
understanding of the stellar content of the cluster and the entire IR-dark filament. Modeling 
of the SEDs of GLIMPSE sources identified 31 young stellar objects in the filament with a 
combined stellar mass of ~ 127±27M . An additional 22 sources were identified as probable 
cluster members based on the presence of strong 24 /zm emission. 

The total star formation efficiency in the G34.4 cloud filament is estimated to be about 
7% while the massive and intermediate mass star formation efficiency in the entire cloud 
filament is estimated to be roughly 2%. 

A qualitative evaluation of the outflow characteristics and presence of ionized gas in 
the central region of the IR-dark cloud filament suggests that the early B stars currently 
undergoing outflow and accretion began forming stars at approximately the same time: 
about 10 5 years ago. A population of lower mass stars were detected by SNB04 that appear 
to have an age of about a few Myrs. This dual population suggests that stars in the IRAS 
18507 cloud may have formed in two stages: first lower mass stars formed and then, a few 
Myrs later, the more massive, early B stars began to form. Such a situation can occur in 
magnetically sub-critical clouds in which slow, ambipolar diffusion leads to the delayed onset 
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of core contraction and eventual massive star formation. 

A comparison of the gravitational binding energy with the outflow kinetic energy in 
the two cloud cores studied in molecular lines suggests that the compact core containing 
G34.4 MM is being destroyed by its molecular outflows. The outflows have as much or more 
kinetic energy than the gravitational binding energy of the cloud core. Thus, we may be 
seeing the first and only episode of massive star formation in this compact region. 

In contrast, the 690 M core surrounding the G34.4 UCHll region has twice the 
gravitational binding energy as the combined outflow kinetic energy. Thus, although the 
outflows are injecting a significant amount of mechanical energy into the cloud core and may 
help prevent further collapse of the cloud, they are not likely to totally disrupt the cloud. 
As long as cloud material remains, this leaves the door open for possible future episodes of 
star formation. 
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Table 1: G34.4 and G34MM Outflow Parameters 



Source: 


G34.4 combined flows 


G34 MM combined flows 


KjVJ IdUlUo Ol UUL11UW 


u.ao pc 


u.oo pc 


.ri.ooU.l.Llt:U. lllUlllldtlUll dJ.J.g.j.t; 






Onfflnw A/TpQet - 
WLilllUW IVIdoo . 






Red-shifted- 12 C0 

j. i v. v.i. oiJ.il w 


22 6 Mr, 


10 6 Mr, 


13 C0 


59.6 M 


12.7 M 


Blue-shifted: 12 C0 


14.2 M 


1.6 M 


13 C0 


14.8 M 


9.9 M 




111.2 M 


34.8 M 


Momentum 


8.3 x 10 2 M km s" 1 


3.5 x 10 2 M km s" 1 


Kinetic Energy 


7.8 x 10 46 ergs 


4.9 x 10 46 ergs 


Dynamical time scale of CO flow 


~ 4 x 10 4 yr 


~ 3.8 x 10 4 yr 


M f 


1.3 x 1(T 3 M yr" 1 


0.8 x 10~ 3 M yr" 1 


Momentum Supply Rate (Force) 


12.2 x 1(T 3 M km s^ 1 yr~ x 


8.0 x 10~ 3 M km s~ x yr" 1 


Mechanical Luminosity 


11.2 L 


9.1 L 



f 12 CO outflow emission (corrected for optical depth from [CO]/[ 13 CO]) 

measured between velocities 32.9 to 51.7kms _1 and 63.6 to 77.8kms~ 1 . 

13 CO outflow emission (assumed to be optically thin) 

measured between velocities 51.7 to 53.6 km s -1 and 60.4 to 63.6kms~ 1 . 



Table 2: G34.4 and G34MM Core Parameters measured from C 18 



Source: 


G34.4 


G34MM 


C 18 Core Mass: 


692.7 M 


75.2 M 


VLSR 


56.0 km s -1 


57.4 km s -1 


Momentum 


2.0 x 10 3 M km s" 1 


2.5 x 10 2 M km s" 1 


Kinetic Energy 


8.0 x 10 46 ergs 


1.2 x 10 46 ergs 
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Figure Captions 

Figure 1. Four-color GLIMPSE image showing 3.6 //m, 4.6 fim, 5.8 /im, and 8.0 /im as blue, 
green, orange, and red colors, respectively. The image is displayed in galactic coordinates 
and is roughly 16.6' x 21.5' in size. The white box at upper left gives the field of view 
shown in Figure 2. The G34.4 UC Hll region is in the center of the box within a dense cloud 
filament seen as a dark lane. The mid-IR counterpart to the G34.4 MM core is seen as a faint 
green smudge immediately to the upper left of the UC Hll region. The larger nebula in the 
lower right of the image is associated with the cometary UC Hll region G34.26. The bright 
green emission emanating from this region traces a very energetic outflow. This GLIMPSE 
image includes the same data as reported in Rathborne et al. (2005) but it covers a wider 
field-of-view to show the physical relationship between the IR-dark cloud in which the G34.4 
UCHll is located and the Hll region G34.26. 

Figure 2. Four images, all showing the same field of view. Clockwise from the upper-left 
panel are: (a) 2MASS K-band, (b) GLIMPSE 4.5 fim, (c) GLIMPSE 8.0 /im, and (d) a 
3-color image combining the three other images as blue, green, and red, respectively. Green 
plus signs mark the positions of the UCHll and G34.4MM, while green circles show the 
locations predicted for the driving sources of outflows G34.4:D and G34.4:E. The white box 
gives the approximate field of view of Figure 3. The locations of the millimeter cores detected 
by Rathborne et al. (2005) in J2000 coordinates are: MM1: 18 h 53 m 18.0 s , 01°25'24"; MM3: 
18 h 53 m 20.4 s , 01°28'23"; MM4: 18 h 53 m 19.0 s , 01 o 24'08"; 

Figure 3. CO(J=l— 0) red and blue-shifted emission. Top left panel shows red and blue- 
shifted emission from 62.3 to 79.2 kms -1 and 37.5 to 53.2 km s" 1 , respectively. The RMS 
is 470 mJy beam -1 kms" 1 in the red-shifted image and 370 mJy beam -1 kms -1 in the blue- 
shifted image. Bottom panels show low-velocity CO emission (left panel, 62.3 to 67.5 kms -1 
and 49.3 to 53.2 kms -1 , RMS = 350 mJy beam -1 kms -1 ) and higher velocity CO emission 
(right panel, 67.5 to 79.2 kms -1 and 37.5 to 49.3 kms -1 , RMS = 350 mJy beam -1 kms -1 ). 
In all contour maps, contours are plotted at 3, 4, 5, 7, 10, 15, 20 a and continue with a spacing 
of 5<T. The top right panel shows the NIR K band image of Shepherd et al. (2004) with 
the locations of the proposed outflows overlaid. In all panels: the solid circle represents the 
location of G34.4 MM, the star shows the location of the UC Hll region, the outflow (A) from 
G34.4MM is outlined with a series of arcs which delineate the approximate boundaries of 
the flow. Four other proposed outflows are shown as arrows. A solid triangle indicates the 
possible location of the driving source of an outflow. The outflow from the UC Hll region 
is one-sided since only the blue-shifted outflow is clearly detected. The synthesized beam of 
3.83" x 3.45" at P.A. -59.7° is shown in the bottom right of the upper left panel. 

Figure 4. C 18 (top, left panel) and 13 CO (top, center panel) integrated emission 
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(moment 0) from 62.1 kms -1 to 54.6 kms -1 and 64.5 km s -1 to 49.5 km s -1 , respectively 
{vlsr — 57kms -1 ). The RMS in the C 18 image is 200mJybeam -1 kms" 1 ; contours 
are plotted at 3, 5, 7, 9, 11, 13, 15, &17<7 while grey scale is shown from 200 mJy beam" 1 
kms -1 to a peak of 3.1 Jy beam -1 kms -1 . The RMS in the 13 CO image is 280 mJy beam -1 
kms -1 ; contours are plotted at 3, 5, 7, 9, 12, 15, 20, 25, & 30<r while grey scale is shown from 
280 mJy beam -1 kms -1 to a peak of 9.0 Jy beam -1 kms -1 . Plus symbols represent the 
locations of the UCHll region and G34.4MM. Green triangles represent the locations of 
spectra taken for optical depth estimates in the outflows. A scale size of 0.55 is shown in 
the left panel, synthesized beams are shown in the bottom right of each panel. The top, 
right panel shows 13 CO red and blue-shifted emission from 63.1 to 61.8 kms -1 and 52.3 to 
49.5 kms -1 , respectively. The RMS in the image is 100 mJy beam -1 kms -1 . Contours are 
plotted at 3, 6, 9, 12, 15, 18 and 21 a. Symbols are the same as in Figure 3. The bottom 
two panels show moment 1 maps of C 18 & 13 CO in which the brightness of the color is 
proportional to the intensity in the moment image. Color bars show velocity and brightness 
scales in each image. 

Figure 5. CO(J=l— 0) channel maps at 1.3 kms -1 spectral resolution. The central velocity 
is indicated in the upper left of each panel. The channel RMS is 60 mJy beam -1 , the peak is 
3.3 Jy beam -1 , and contours are plotted at ±5, 15, 20, 30, 40 and 50 a. The lower right panel 
shows the synthesized beam in the bottom right corner (3.83" x 3.45" at P.A. —59.7°) and a 
scale size of 0.55 pc. Plus symbols represent the locations of the UCHll region at the 6 cm 
continuum peak and G34.4 MM as indicated in the upper left panel. 

Figure 6. 13 CO(J=l— 0) channel maps at 1.36 kms -1 spectral resolution. The central 
velocity is indicated in the upper left of each panel. The channel RMS is 44mJybeam -1 , 
the peak is 1.7 Jy beam -1 , and contours are plotted at ±3,5,7,9,12,15,20 and 25 a. The 
lower right panel shows the synthesized beam in the bottom right corner (4.04" x 3.83" at 
P.A. —49.9°) and a scale size of 0.55 pc. Plus symbols represent the locations of the UCHll 
region at the 6 cm continuum peak and G34.4 MM as indicated in the upper left panel. 

Figure 7. C 18 0(J=1— 0) channel maps at 0.7 kms -1 spectral resolution. The central 
velocity is indicated in the upper left of each panel. The channel RMS is 60mJybeam -1 , 
the peak is 0.92 Jy beam -1 , and contours are plotted at ±3, 5, 7, 9, 11, 13 and 15 a. The lower 
right panel shows the synthesized beam in the bottom right corner (4.04" x 3.83" at P.A. 
—50.1°) and a scale size of 0.55 pc. Plus symbols represent the locations of the UC Hll region 
at the 6 cm continuum peak and G34.4 MM as indicated in the upper left panel. 

Figure 8. Spectra convolved with a 10" beam taken in CO, 13 CO and C 18 at different 
locations in the image. The ratio of the CO and 13 CO emission was used to calculate optical 
depth in the red and blue-shifted outflow lobes as indicated by the dashed-dot lines and 



- 32 - 



arrows in the figure. Panels i) & ii) show spectra taken near the G34.4 UCHll region at 
positions 18 h 53 m 18.320 s 1°24'45.84" and 18 h 53 m 19.203 s 1 24'39.38", respectively. Panels iii) 
& iv) show spectra taken near the G34MM core at positions 18 h 53 m 18.065 s l o 25'25.06" and 
18 h 53 m 17.284 s 1 25'13.52", respectively. The locations of the spectra relative to the 13 CO 
emission is shown in the center panel of Figure 4. The dashed vertical line at 57kms _1 helps 
to illustrate that the peak of the C 18 line shifts from 57.4 km s -1 toward the G34MM core 
to 56kms _1 toward the G34.4 UCHll region core. 

Figure 9. GLIMPSE 4.5 /zm image in inverted greyscale showing locations of G34.4MM 
and G34.4 UCHll region as crosses. Open diamonds show the predicted locations of the 
driving sources of outflows G34.4:D and E based on CO and 13 CO morphology. GLIMPSE 
sources are located very near the predicted positions suggesting that GLIMPSE has detected 
two sources driving outflows D and E. 13 CO red- and blue-shifted emission contours from 
Fig. 4 are overlaid for reference. 

Figure 10. GLIMPSE and MIPSGAL 3-color image with 4.5 ^m (blue), 8.0 fjm (green) and 
24 /im (red). The probable YSO cluster members listed in Table 3 are marked by color-coded 
circles: green = good YSO fits (x 2 /N data < 4.8), yellow = bad YSO fits (x 2 /N data > 14.8), 
cyan = source detected in fewer than four bands, and green = MIPS 24yum point sources 
detected in too few IRAC+2MASS bands to be fit. Sources well-fit by stellar photospheres 
are marked by blue circles. The IR dark cloud has been divided into 3 regions for the 
purpose of the YSO census: Region C (large white circle) contains the OVRO observations 
and regions N and S are (large white rectangles) are the portions of the dark cloud lying to 
the north and south, respectively, of region C. Source numbers from Table 3 for the YSOs 
discussed in Section 3.4 are also shown. 

Figure 11: (a) Ten best YSO SED model fits to the IR fluxes of source 12/G34.4MM. The 
fitting assumes that the mid-IR flux from the millimeter core is dominated by a single source. 
The [4.5] upper limit (inverted triangle) was used because the source appears to have excess 
flux in this bandpass due to molecular emission associated with outflow activity that has not 
been incorporated into the YSO models. The thickest line represents the best fit to the data, 
and different colors represent the SED observed through different apertures, from smallest 
(IRAC, black) to largest (MSX 21.3 fim, green). The dashed line shows the photospheric 
flux from the best-fit YSO model, as it would appear if suffering from interstellar extinction 
only in the absence of circumstellar material, (b) Masses and luminosities of the best fit 
models. The shaded contour regions represent the logarithmic density of models in the 
mass-luminosity parameter space. The darkest gray represents a model density 1000 times 
higher than that of the lightest grey. The white region beyond the lowest contour shows the 
part of parameter space that lies beyond the edge of the model grid. 
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Figure 12: (a) Same as Figure 11a, but for the other 6 well-fit YSOs in region C. All fits 
used to compute the model parameter ranges in Table 4 are plotted. Sources 10 and 11 
overlap the UC Hll region, and hence the diffuse mid-IR fluxes associated with the UC Hll 
region in all 4 MSX bands and MIPS 24 and 70 /im were used as upper limits to the fits, 
(b) Same as Figure lib, but for the 6 well- fit YSOs plotted in Figure 12a. 
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Figure 6. 
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Figure 7. 
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Figure 8. 
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Figure 9. 
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Figure 10. 
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Figure 11a. 
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Figure 12a. 
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Figure 12b. 



